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Both the physical picture of the dynamics of atoms and molecules in intense infrared fields and its 
theoretical description use the concept of electron trajectories. Here we address a key question which 
arises in this context: Are distinctly quantum features of these trajectories, such as the complex-valued 
coordinates, physically relevant in the classically allowed region of phase space, and what is their origin? 
First, we argue that solutions of classical equations of motion can account for quantum effects. To this end, 
we construct an exact solution to the classical Hamilton-Jacobi equation which accounts for dynamics of 
the wave packet, and show that this solution is physically correct in the limit h —> 0. Second, we show 
that imaginary components of classical trajectories are directly linked to the finite size of the initial 
wavepacket in momentum space. This way, if the electronic wavepacket produced by optical tunneling in 
strong infrared hels is localised both in coordinate and momentum, its motion after tunneling ipso facto 
cannot be described with purely classical trajectories - in contrast to popular models in the literature. 


1. Introduction 


Nonclassical solutions to the Hamilton-Jacobi equation (HJE) [lj, known as quantum trajectories, 
play an important role in strong-field physics 0-0. They allow one to seamlessly incorporate 
quantum aspects of the dynamics, notably the laser-induced tunnelling, keeping at the same time 
a physical insight into the problem. 

The assumption underlying many theoretical models of strong-field ionization dynamics - from 
the famous three-step model 0,0 to Coulomb-corrected models of strong-field ionization Ml 
to the analysis of attoclock experiments fTTilfil] - is that quantum trajectories are limited to the 
motion in imaginary time in the classically forbidden region during the tunnelling stage. So far 
as the electron has emerged from under the barrier, its motion to the detector is expected to be 
classical in the true sense of the term. Consequently, in simulations, both the coordinate and the 
momentum characteristic of its trajectory are presumed to be real [l2-15|]. 

The notable exception are the studies of high harmonic generation and correlated double ionisa¬ 
tion, where the electron momentum remains complex after ionisation, until its recollision with the 
parent ion (see, e.g., [l6[ and references therein). Furthermore, it has recently been shown 17-21] 
that the presence of an imaginary component in an otherwise classically-behaving trajectory af¬ 
ter tunnelling is crucial for quantitative agreement with ab-initio calculations even in the absense 
of recollision (2^|. While for short-range potentials imaginary part of the trajectory coordinate 
remains hidden from view, as it does not impact ionization yields or photoelectron spectra, for 
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long-range potentials it leads to qualitative and quantitative consequences f23j. Counter-intuitive 
as they may seem, quantum trajectories appear to exist in any case of classically permitted motion, 
be it from tunneling to recollision or from tunneling to the detector. 

We argue that resolution of this controversy is in the fact that solutions to classical equations of 
motion, such as the HJE, need not be by themselves classical. They can account to some extent 
for quantum effects. Emergence of quantum trajectories in classically permitted regions is thus 
associated with a certain hierarchy of quantum effects, with some accounted for in solutions to the 
HJE, and others neglected. 

In this note, we demonstrate real-time quantum trajectories for the simplest textbook exam¬ 
ple: the motion of a Gaussian wavepacket in a field of homogeneous time-dependent force. The 
assumption of a Gaussian wavepacket is typical, e.g., in the analyses of the electron dynamics 
after tunneling induced by the infrared laser pulses. Following the standard logic of separation of 
“exponential terms” from “pre-exponential factors”, we regard normalisation of the wavepacket 
“pre-exponential”. We show that this is justified in the semiclassical limit H —> 0. The “exponential 
terms” then add up to an exact nonclassical solution to the HJE. Quantum interference (phase) 
effects, which are, according to such logic, “exponential”, are by construction included in this so¬ 
lution. Nonclassical nature (complexity) of the latter is shown to be associated with localisation of 
the Gaussian packet in coordinate and momentum. 


2. WKB Ansatz revisited 


2.1. WKB basics 

Formally, the aforementioned exact nonclassical solution to the HJE emerges within a quasiclassical 
(WKB) approximation 24-26j to evolution of a Gaussian wave packet under the influence of a 
homogeneous force. Under the said approximation, one solves the Schrodinger equation, 


ihd t ^{r,t) = 


~ 77-dr +V( r it) 
2m 


(1) 


by making the Ansatz, 


t ) = exp 


iS(r, t ) 


+ ia(r,t ) 


The “WKB action” S(r,t ) obeys the equation, 

—d t S(r , t) = d r S(r,t )] 2 + V(r,t), 

which is a classical HJE for the Hamiltonian, 

H(p,r,t) = %— + V(r,t). 

2m 

Substituting this to (P) we find the equation for a(r,t), 


( 2 ) 

( 3 ) 

( 4 ) 


- d t a(r,t ) 


— d r S(r , t)d r a(r,t) 
m . 


i AS(r '* ) + 2^ {|9 ’" T(r ' t)l! 


iAa(r, f)}. 


( 5 ) 


In the quasiclassical (WKB) approximation, the term proportional to H in flSJ) is dropped. 
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2.2. Is WKB action necessarily classical? 

Equations (JH]) , © follow by substituting the WKB Ansatz © in the Schrodinger equation (pQ) and 
grouping the terms according to powers of h. However, this is no more than a leading consideration. 
By themselves, Eqs. ©, © are exact. More precizely speaking, with the HJE © postulated, 
Eqs. m and © are equivalent. 

There is nothing formally wrong with accommodating quantum effects in S(r,t). Such approach 
may be beneficial if a particular quantum effect exceeds in magnitude certain classical ones. Indeed, 
think of a wave packet propagating in a potential. Such quantum motion can always be separated 
into the classical (trajectorial) dynamics and evolution of the wave-packet. If the potential is 
linear (a homogeneous force), the packet dynamics reduces to spread. For a potential that is 
only approximately linear (a nearly homogeneous force), the quantum spread may well happen to 
remain large compared to corrections to the classical trajectory. In this case, classification of effects 
according to their classical or quantum nature does not make much sense. 


2.3. WKB action versus the classical action 

Now, where is the mathematical opening for S(r,t ) to include quantum effects? It is instructive to 
compare S(r,t ) to the classical action proper, 


S(r,t-r 0 ,t 0 )= f L(p(t'),r(t'),t')dt', (6) 

Jto 

where the classical trajectory p(t),r(t) is specified by the conditions, 

r(t) = r , r(t 0 ) = r 0 . (7) 

The classical action is a “two-ended” quantity. It is specified by two boundary conditions for 
trajectories, and obeys in fact two HJE’s, Eq. © and the other one with dt —> —dt 0 , d r —> —<9 ro , 
V(r,t ) —> K(ro,to)- The WKB action S(r,t ) is a “one-ended” quantity. It is specified by one 
boundary condition , r(t) = r, and obeys one HJE ©. That is, conditions specifying the WKB 
action are weaker than those for the classical action in the true meaning of the term. 

To see how this opening may work, consider restoration of quantum trajectories from a given 
WKB action. They are found from the equation, 

md t r=p(r,t), ( 8 ) 

where p(r,t) is the momentum as a function of r,t, 

P{r,t) = d r S(r,t). (9) 

Equation © is a system of three first-order differential equations; their solution depends on three 
integration constants. This fact is of utter importance. It shows that by specifying action one 
reduces the number of available integration constants from six - cf. Eq. © - to three, thus impos¬ 
ing a condition on underlying trajectories. For the classical action, this is the boundary condition 
r(to) = tq. However, in general, there is no guarantie that this condition may be interpreted clas¬ 
sically. In turn, this would mean that the WKB action is itself a nonclassical quantity - the fact 
that it obeys the HJE notwithstanding. 
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3. Nonclassical solutions to the Hamilton-Jacobi equation 
3.1. Gaussian wave packet in a field of homogeneous force 

A physically meaningful example of a nonclassical WKB action emerges if considering evolution of 
a Gaussian wave packet in an arbitrary nonstationary homogeneous force field, 

V(r,t) = V L (r,t) = -rF(t). (10) 

An exact solution to this problem reads, 

3/4 exp (1 {wl W 4- [r - X(0]r(t)} - (11) 

Here 

D(t) = D(t 0 ) + (12) 

m 

determines the width of the wave packet. Other quantities in (Hill are defined for a classical motion 
in potential (flOl) : X(t),P(t) is a classical trajectory obeying the Hamilton equations, 

d t x(t) = d t P(t) = F(t), (13) 

m 

and 

P P' 2 (t') 

W L (t) = X(t)P(t) - jf (14) 

The intial conditions X(to), P(t o) and D(to) may be arbitrary, except the real part of D(to) must 
be positive to assure normalisation of (HU), 

9£D(t 0 ) > 0. (15) 

That (fITIl solves the Schrodinger equation (P) with V = Vl is readily verified by direct calculation. 


f(r,t) = 


&D(t 0 ) 

irD 2 (t,) 


3.2. Nonclassical solution to the HJE 

Solution dill) should be compared to the WKB Ansatz Q. The obvious choice of the WKB action 
S(r,t ) is to identify it with the combination in curly brackets. A less obvious choice is to account 
in S(r,t ) for the quantum phase stucture this packet acquires when spreading, 

ih\v — 

S(r,t ) = W L (t) + [r - X(t)]P(t) + [ 2D{t) [ . (16) 

The “quantum correction” is then limited to normalisation, 

l o(r,t)=io(t) = -ln-p^. (17) 
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To verify that ([TO indeed obeys the HJE, 

dtS{r,t) + -!-[ d r S(r,t )] 2 - rF{t ) = 0, 


(18) 


is again straightforward. 


3.3. Quantum trajectories 

So what kind of motion does S(r,t ) describe? The momentum as a function of r,t reads, 


P(r,t) = P(t ) + 


ih[r — X(t)] 

W) 


Solving equation @) we obtain the trajectories, 


where 


r(f) = X{t) + 


D(t)Pi 

ih 


pit) = P(t) + Pi, 


ih[r{t 0 ) - Xq(Q| 

D(to ) 


(19) 


( 20 ) 


( 21 ) 


is the triad of integration constants we spoke about in Sec. 12.31 
While functions (1201) do obey classical equations of motion, they are not what one would normally 
call “a classical trajectory”. To start with, they are complex. Reality of the momentum may be 
maintained (with real p,), but not reality of the coordinate. These peculiarities are rooted in the 
initial condition, 


ih[r(t 0 ) - X(t 0 )] = D(t 0 )\p(t 0 ) - P(t 0 )]. (22) 

It makes no sense in classical mechanics because of restriction (|15l) . We stress that time t remains 
real, so that the WKB action S(r,t) and the trajectory (EU1) are characteristic of the classically 
permitted motion. 


4. Discussion 

The WKB action (fTHI) emerged as an ad hoc quantity. In fact its choice may be supported by 
scaling arguments in the classical limit h —> 0. Assume that D(t) is real at some time t i n . For large 
propagation times, the wave packet spread becomes much larger than its initial width (cf. Sec. \2.2\i . 
Neglecting the latter, 


and 


D(t) 


ih(t - tin) 

m 


S{r , t) 


W L (t) + [r 


X(t)]P(t) + 


m[r — X(t)] 2 
2 (t - tin) 


(23) 


(24) 
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Planck’s constant has dropped out. In the approximation ()23l) . it only survives in the preexponential 
normalisation factor. Influence of this factor in the limit h —> 0 is small. This is clear from comparing 
S(r,t)/h to a(r,t). The former scales as 1/h while the latter as Info. 

It is interesting and instructive to note that in the approximation (1231) trajectories (1201) and the 
condition (122[) become real, and start making formal sense in classical mechanics (with real pf). In 
particular, Eq. (l22|) becomes, 


r(t 0 ) - X(t 0 ) = ——— [p(t 0 ) - P(to)]. 


m 


This condition makes perfect classical sense. Recalling Eq. (1201) it may be rewritten as, 

(to ~ tin)Pi 


r(t 0 ) - X(t Q ) = 


m 


(25) 


(26) 


or simply 


r(t in ) = X(t- m ). (27) 

That is, in the approximation (l23l) . WKB action (flfil) reduces to the classical action (j6j), with 
to —> ti n and tq —>• X(ti n ). Up to a constant, 

S(r,t) = S(r,t-,X(tin),tin). (28) 


Physically, it describes a bunch of trajectories “growing” from a single point X(t- m ), quite natural 
a result of neglecting the initial width of the packet. 

This way, nonclassically of the action (1161) hinges on whether the initial width of the Gaussian 
pachet can or cannot be neglected. In any case, quantum mechanics prohibits Eq. (1231) as an exact 
relation, because it would render the wave packet unnormalisable. Complexity in (1201) and (|22l) are 
thus traces of quantum mechanics in “nearly classical” formulae. 

This is the case, for example, in interpretation of the attoclock experiments on tunneling. The 
deviation of the observed results from those based on classical modeling of the electron dynamics 
after ionization are attributed to the parameters of the initial Gaussian wavepacket assumed after 
tunneling. Once these parameters become important, the quantumness of the underlying trajecto¬ 
ries has to be included. The quantitative consequences of this quantumness have been discussed in 
0 - 


In conclusion, a few words so as to put this paper in a wider perspective. After the initial shock 
- how on earth the electron choses the state to leap to - there has been a growing awareness that 
quantum and classical mechanics overlap rather than are mutually exclusive. Informally speaking, 
the quantum-classical border is an area rather than a line. This viewpoint has been especially 
fruitful in quantum optics, culminating in Glauber’s photodetection theory 27-0 and phase- 
space techniques (3CH37l|. 

One of the defining features of the said overlap is that quantum processes may be governed 
by classical equations of motion. For the electroma gne tic field, this extends in fact to full formal 
identity of the quantum and classical dynamics 38j, l39|- Ipso facto , such identity is impossible for 
equations of motion of matter, because it would mean rewriting the whole of quantum mechanics 
as a hidden-variable theory |40|] (and this is not to mention fermions). The critical question is 
therefore, how far classical concepts actually extend into the realm of quantum mechanics. As 
is shown in this paper, one of the “inhabitants” of the quantum-classical overlap is the classical 
Hamilton-Jacobi equation. 

To summarise, the classical action is a solution to the HJE, but not vice versa. Even in real time, 
an arbitrary solution to the HJE needs not be classical. Specifying it implicitly imposes conditions 
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on “classical trajectories”, which may happen to be impossible for classical trajectories in the true 
meaning of the term. One may say that the Hamilton-Jacobi equation is in fact a quantum equation, 
of which some (but not all) solutions make sense in classical mechanics. “Quantum” solutions to 
the HJE are physically important within the quantum-classical overlap, i.e., in the limit h —> 0. 
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